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We use a ferromagnetic voltage probe model to study the influence of inelastic scattering on giant 
magnetoresistance and current-induced torques in ferromagnetic and antiferromagnetic metal spin 
valves. The model is based on the Green's function formulation of transport theory and repre- 
sents spin-dependent and spin-independent inelastic scatterers by interior voltage probes that are 
constrained to carry respectively no charge current and no spin or charge current. We find that 
giant magnetoresistance and spin transfer torques in ferromagnetic metal spin valve structures sur- 
vive arbitrarily strong spin-independent inelastic scattering, while the recently predicted analogous 
phenomena in antiferromagnetic metal spin valves are partially suppressed. We use toy-model nu- 
merical calculations to estimate spacer layer thickness requirements for room temperature operation 
of antiferromagnetic metal spin valves. 

PACS numbers: 



I. INTRODUCTION 

Electronic phase coherence often plays an important 
role in mesoscopic quantum transport. The simplest ex- 
ample, perhaps, is the Aharonov-Bohm effect in a meso- 
scopic ring which is manifested by dependence of con- 
ductance on enclosed magnetic flux and requires phase 
coherence across the ring. In the Landauer-Biittiker the- 
ory of quantum transports 2 * 3 - it is possible to simulate 
phase breaking scattering processes by including inter- 
nal 'floating' voltage probes in addition to source and 
drain^ electrodes, and requiring that their chemical po- 
tentials adjust so that they do not carry a current. This 
simple model is compatible with a Green's function de- 
scription of transport and is able to describe those qual- 
itative consequences of inelastic scattering not associ- 
ated with the peculiarities of a specific mechanism. It 
has, for example, been applied to study the influence 
of phase decoherence on spin-independent transport in 
one-dimensional systems*^, and the effect of dephasing 
on transport through quantum dotsi*V 

In this paper we report on a comparison of the influ- 
ence of inelastic scattering on giant magnetoresistance 
(GMR) and spin-transfer (ST) torque effects in tra- 
ditional spin-valve structures containing ferromagnetic 
metals, and in the antiferromagnetic spin-valve struc- 
tures that w o?' 1 ? have recently proposed. (Spin- valves 
are illustrated schematically in Fig. QJ) These devices 
are normally intended for operation at room tempera- 
ture and above and electron scattering is therefore nor- 
mally dominantly inelastic. In traditional ferromag- 
netic spin-valves only spin-flip scattering in the para- 
magnetic spacer is expected to have a large impact on 
GMRpiSJi or STi-U£. Both effects are expected to be 
strong when the spin diffusion length is longer than the 
paramagnetic spacer thickness. (Spin-independent in- 



elastic scattering does however play a role in limiting 
the amplitude of oscillations in the dependence of ex- 
change coupling between ferromagnetic layers on spacer 
layer thickness^.) It is commonly believed that GMR 
and ST in ferromagnetic spin-valves can survive arbi- 
trarily strong spin-independent inelastic scattering. For 
antiferromagnetic spin- valves, however, the analogous ef- 
fects depend^, at least in a simple toy model, on coher- 
ent multiple-scattering in the spacer layer and at its in- 
terfaces with the antiferromagnets. We have therefore 
argued that both effects will become weak when inelas- 
tic scattering is strong, even if the scattering is spin- 
independent. Since phonon scattering is only weakly 
spin-dependent and strong at high temperatures, inelas- 
tic scattering cannot be ignored in practical antiferro- 
magnetic spin valves. 

To pursue these ideas more quantitatively we use a 
voltage probe model to represent the influence of spin- 
independent and spin-dependent inelastic scatterers in a 
Green's function description of transport through a spin 
valve. In the linear response regime the voltage probe 
model has been shown to be equivalent to local cou- 
pling to a harmonic-oscillator heat bathii. In a mag- 
netic metal circuit, however, conventional voltage probe 
models produce both longitudinal and transverse spin re- 
laxation. Consider for example a largely collinear mag- 
netic configuration with a natural spin-quantization axis. 
Spin components transverse to this axis are represented 
in quantum mechanics by coherence between spin-up and 
spin-down projections of the form | f) +e l ^ \ |). A conven- 
tional paramagnetic voltage probe will change the mag- 
nitude and randomize the phase of both | f ) and | |) 
parts of this spinor independently and therefore alter all 
spin-density components. This property was used very 
recently by Michaelis and Beenakker to model spin de- 
cay in chaotic quantum dotsii. This feature is however 
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FIG. 1: a) Parallel configuration and equivalent circuit consisting of two resistors R and r < R, and b) antiparallel configuration 
and equivalent circuit, of a spin valve consisting of two single-domain ferromagnets separated by a paramagnetic spacer, 
c) Antiparallel configuration and equivalent circuit of an antiferromagnetic spin valve and d) its parallel configuration and 
equivalent circuit. (For convenience, only two ferromagnetic layers within each antiferromagnet are shown.) Note that in 
the ferromagnetic case the resistance of parallel and antiparallel circuits is different, contrary to the equivalent circuits in the 
antiferromagnetic case. 



undesirable in modelling the influence of inelastic scat- 
tering on magnetoresistance properties of spintronic de- 
vices, because the strongest inelastic scatterers are of- 
ten phonons, and phonons conserve spin to a good ap- 
proximation. In our calculations we therefore use a volt- 
age probe model generalized to the case of ferromagnetic 
probes, shown in Fig. [21 This generalization allows us to 



consider separately spin-dependent and spin-independent 
inelastic scatterers. The spin-independent case is realized 
by adjusting the magnetization direction and the major- 
ity and minority spin chemical potentials of the voltage 
probes so that the probe carries neither charge or spin 
currents. 



There are a number of important distinctions be- 
tween GMR and ST physics in ferromagnetic and 
antiferromagnetic? structures. In the ferromagnetic case 
GMR can be understood qualitatively simply by using 
a two-channel resistor model illustrated in Fig. ^ which 
does not rely on phase coherence. When applied to an- 
tiferromagnetic spin valves, the same argument does not 
predict GMR because the two spin channels have identi- 
cal conductances (see Fig. ^ Clearly the antiferromag- 
netic spin-valve GMR effect relies at least in part on 
phase-coherence near the spacer layer as. In the sim- 
ple model we discuss below antiferromagnetic spin-valve 
GMR comes about, roughly speaking, because the re- 
flection amplitude for an electron scattering off an an- 
tiferromagnet is spin dependent, although the reflection 
probability and transmission amplitude and probability 
are not. Similarly ferromagnetic spin- valve ST torques 
can be understood by using a model of a ferromagnetic 
metal which has quasiparticle and magnetization orienta- 
tion degrees of freedom, and appealing to conservation of 
total spin angular momentum to infer an action-reaction 
relationship between the torque exerted on the quasipar- 
ticles by the exchange field and the torque exerted on the 
magnetization by current-carrying quasiparticles. In the 
antiferromagnetic case, it is the staggered moment anti- 
ferromagnetic order which behaves collectively. Since this 



coordinate does not carry total spin, its current-driven 
dynamics is not specified by a global conservation law. 
The theory of current-driven order parameter dynamics 
in this case requires a more microscopic approach^ in 
which ST is seen as following from changes in the spin- 
dependent exchange potential experienced by all quasi- 
particles, which follow in turn from changes in the spin- 
density in the presence of a transport current. Indeed 
the use of the term spin transfer torque is perhaps in- 
appropriate in the antiferromagnetic case since current- 
induced order parameter changes are not related to trans- 
fer of total spin angular momentum between subsystems. 
Nevertheless staggered torques do act on the staggered 
spins and do drive drive the antiferromagnetic order pa- 
rameter. To make clear that current-induced order pa- 
rameter changes in antiferromagnets are not related to 
conservation of spin and are therefore strictly speaking 
not an example of spin transfer we, from now on, will 
call these torques current-induced torques (CIT's). Re- 
markably, current-induced torques in an antiferromag- 
netic spin- valve act throughout the volume of the antifer- 
romagnets in the absence of inelastic scattering and dis- 
order. Critical currents for order parameter reversal are 
therefore independent of antiferromagnetic film thickness 
in this limit, rather than being inversely proportional to 
thickness as in the ferromagnetic case. 




FIG. 2: Model system consisting of sites with an arbitrary onsite exchange potential. In addition to the coupling to its nearest 
neighbors, each site is coupled to a ferromagnetic reservoir. 



The main aim of this paper is to shed light on the 
robustness of these GMR and current-induced torques 
in realistic room-temperature thin-film antiferromag- 
netic spin-valve structures by examining how properties 
change upon introduction of spin-dependent and spin- 
independent inelastic scatterers. The ferromagnetic volt- 
age probe models of elastic and inelastic scattering are 
completely satisfactory for this purpose because we are 
interested in achieving a qualitative understanding that 
transcends the details of specific systems. We solve the 
model using the Green's function formalism for electronic 
transport in mesoscopic systemsi^iSi. This is convenient 
because the spin-torques we want to evaluate can be 
expressed^ in terms of the transport steady-state elec- 
tron spin density and local observables are readily calcu- 
lated using the Green's function formalism. 

Our paper is organized as follows. In Sec. [H] we 
describe the Green's function formalism for electronic 



transport as applied to systems with ferroma gnet ic leads. 
(For related work see Refs. and Ref. The 
main result of this section is a general expression for the 
spin current from the ferromagnetic leads into the system 
which is necessary to apply our voltage probe model of a 
spin-independent inelastic scatterer. In Sec. lIIII we apply 
the formalism to study the effect of spin-independent and 
spin-dependent inelastic scattering on magnetoresistancc 
and current-induced torques in both ferromagnetic and 
antifcrromagnetic systems. We conclude that although 
GMR and current-induced torques in antiferromagnetic 
spin-valves will be weakened by the inelastic scattering 
always present at spintronic device operation tempera- 
tures, the effects predicted in Ref. should still be eas- 
ily observable in favorable materials. We end in Sec. IIVI 
with a discussion of our results and some suggestions for 
future work. 



II. GREEN'S FUNCTION FORMALISM WITH 
FERROMAGNETIC LEADS 

In this section we describe our model system and give 
an expression for the spin current from the leads into the 
system. 



A. Model 

We consider first a one-dimensional non-interacting 
electron system that is connected on every site j to a fer- 



romagnetic reservoir at equilibrium with chemical poten- 
tial fJ (see Fig. 01 . The non-interacting electrons should 
be understood as quasiparticles in a mean-field descrip- 
tion of a magnetic metal similar to the Kohn-Sham quasi- 
particles of spin-density functional theory. The assump- 
tion of one-dimensionality is not essential. Calculations 
similar to the ones we describe which allow for a number 
of transverse channels could be carried out to model par- 
ticular materials systems realistically. The generalization 
of the formalism to two and three dimensions is straight- 
forward. The total Hamiltonian TL — Hs + "Hl + 7~Lc is 
the sum of three parts. The first term describes the sys- 
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tern, which in our case includes the paramagnetic spacer 
layer and the magnets. In terms of the second-quantized 
operators c J)CT that annihilate an electron in spin state 
a e {T, 1} at site j 



Ti-s = -Js E £ ia £ i'>° 

+ E £ L ( v j s <w - A " n f • t ^-) ^ 



(1) 



In this expression the first term is proportional to the 
hopping amplitude Js and its sum is over nearest neigh- 
bors only. In the second and third terms of 7is we al- 
low for a site-dependent scalar potential Vf and a site- 
dependent exchange potential — A^n| where is a 
unit vector which specifies the (instantaneous) magne- 
tization orientation on site j and t is the Pauli spin- 
matrix vector. This single-particle Hamiltonian should 
be understood as a time-dependent mean-field Hamilto- 
nian that depends on the instantaneous fl^ values. We 
assume that the magnetization dynamics is always slow 
enough to justify time-independent quasiparticle trans- 
port theory. The leads are described by the Hamiltonian 
Hl = X) j 7~t Lj ) where T~t Lj is the Hamiltonian of the lead 
that is connected to the system at site j, which is given 
by 



^ - E [fa 

V 



a j",<? 



r,c>fa , (2) 



where dj 3 a are the fermionic annihilation operators of the 
j-th lead and J-^ j is the hopping amplitude in that lead. 
Each lead has a chemical potential \ij and a uniform ex- 
change potential — A Lj fi Li which we adjust as described 
earlier to simulate spin-independent and spin-dependent 
inelastic scatterers. Note that f2 Lj is in general differ- 
ent for different leads. The Hamiltonian that couples the 
system and its leads is 7ic = Ylj where 



= -4E 



J j,<J U dL j ,a 



(3) 



Here, dLj denotes the last site of the half-infinite fer- 
romagnetic reservoir connected to site j, and Jq is the 
amplitude to hop from that reservoir to the system. 



B. Quantum Transport Green's Function 
Formalism 

The quantum transport Green's function formalism 
determines the equal-time "lesser" Green's functioni^^ 



(t)cj t<T (t)) of the system in the 



transport steady state, the quantity in terms of which 
all observables are calculated. It is given by 1 -' 2 - 



)Ai(e) , (4) 



where N(x) = {eP x + 1 ) 1 is the Fermi distribution func- 
tion, (3—1/ k#T the inverse thermal energy, and 



A j (e) = G^(e)hT j (e)G^(e) , 



(5) 



the spectral weight contribution from lead j. The matrix 
elements of the rate HT^(e) are related to the retarded 
self energy of lead j, denoted by fi,E J, ( + ) (e), by 



tiT j (e) = i \hY, j ' i+ \e) - h^^(e) 



(6) 



where HYP'^^e) is the complex transpose of the retarded 
self energy. (G^- 1 is obtained similarly from the expres- 
sion for G(+) below.) The only nonzero elements of the 
retarded self-energy are 



H& (+ ) ,(e) 



ikl (e)a 



2J Lj 



2Jl 

ikl (e)a ikl (e)a 

e t w — e i w 



8<j,cr> 



(7) 



where fcj(e) = arccos[— (e + aA Lj •)/2Jj J ]/a is the wave 
vector in the leads at energy e, with a the lattice constant. 
Finally, the retarded Green's function is specified by 



G(+)(e) = l, (8) 



with e + — e + iO where the system Hamiltonian is 



H-i 



-WOf'W A^r-T,^,, ) . (9) 



Note that H, p, G<, A j , hT j , and G^ are ma- 

trices in real and spin space of dimension (2iV s ) x (2N S ), 
where N s is the number of sites in the system. From now 
on we use the convention that any quantity is a matrix 
in those indices that are not explicitly indicated. For 
example, the quantity (e) is a 2 x 2 matrix in spin 
space. 

We now proceed to calculate the current and spin cur- 
rents from the leads into the system. Since the underly- 
ing model Hamiltonian in these model calculations are 
spin-rotationally invariant, each component of spin is 
conserved and there is no ambiguity in the spin-current 
definition. The expression for the current from lead j has 
been derived previously^ and is given by 



dW 
~dT 



t 



.(*) 
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Tr 



nP(e)GW(e)nr fc (e)G ( " ) (e)l ■ (10) 



This expression relates the Landauer-Biitikker formalism 
to the Green's function formalism^!. 

The expression for the spin current into lead j proceeds 
along similar lines. We first calculate the rate of change 
of spin-density in lead j: 



~~dt ~ dt 



tz %^(*) 



2h 



j';<r,<r' 



.(ii) 



In principle, there is an additional term that describes 
the precession of electronic spins in the lead around the 
exchange field in the leads, but this term vanishes since 
the leads are assummed to be in equilibrium and assumed 
to have a spin that is aligned with the exchange field. 

The evaluation of the above expression follows the 
same line as in the charge current expression derived by 
Meir and WingreenS 3 .. We first introduce the Green's 



function G^ a ,(t,t') = i 



The 



corresponding Keldysh contour ordered Green's function 
obeys the Dyson equation 23 

Gf(t,t') = £f dt"G jd (t,t")G^ dh .(t",t') , (12) 

where the time integration is over the Keldysh contour 
C°°. In Eq. (|T2"j) G(t,t') is the contour ordered Green's 



function of the system coupled to the leads whereas, in 
order to avoid double-counting of the effects of coupling 
between system and leads, G Lj (t, t') is the contour or- 
dered Green's functions for the j-th lead in the absence 
of coupling to the system. Using standard rules for cal- 
culus on the Keldysh contour we find that 



Gf <<(*,*') - ^ I dt" G^>(t,t")G^ dL] (t",t') 



V 



(+)/ 



+ G< 3 (t,t")G^ 3 {t",t') 



which after Fourier transformation results in 



Gf<{e) = Jl G^(e)G L £< dLj (e) 



+Gfj(e)G h £^ L (e) 



This expression is evaluated using 



G 



L 3 -,< 



(e) 



GdL^dLj ( £ ) 



(4) 



N{e-^)hT] tj {e) 



(e) 



and the kinetic equation [Eq. J3J] to obtain 



X;^(e-At*)4j(e)fiB 



j.(-) 

j,3 



It follows that 



(13) 



(14) 

(15) 
(16) 



(17) 



dsi 
~~dt 



2h 



dt 
(2n 



rG^ie) 



G { ~ht)T 



V 



k 



N(e~^i k )i 



the 



(e) 



h^f\e)T 



A 



In the next subsection we turn to a discussion of the 
physical content of this equation. 



C. Spin-Currents, Exchange Interactions, and 
Current-Induced Torques 

In order to provide some insight into the general ex- 
pression we have derived for the spin-current, we examine 
first the relatively simple situation in which the system 
is connected to leads only on the most left and most 
right sites, denoted by site 1 and N s respectively. Their 
respective chemical potentials are \i x — ep + |e|V, and 



H Ns = ep with \e\V > 0. Moreover, we assume that there 
is no exchange potential in the system, i.e., A| = 0, 
so that spin currents are conserved in the system. This 
is the circumstance then of transport between ferromag- 
netic leads through a paramagnetic system. In this case 
ds 1 /dt = —ds N "/dt. It is informative to separate the 
spin-current flowing between leads into equilibrium and 
non-equilibrium contributions: 



~dt 



ds 1 
~~dT 



cq 



ds 1 
~dt 



(19) 



acq 



It follows quite generally from the time-dependent 
Schroedinger equation satisified by quasiparticles that 
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ds/dt has contributions from the spin-current divergence 
and from precession around an effective magnetic field. 
Since the leads are by definition spatially homogeneous, it 
follows that spin-currents in the leads may be equally well 



thought of as a torque acting on the leads. The equilib- 
rium spin-torque is given (at temperatures low compared 
to Fermi energies) by 



ds 1 
~~dt 



cq 



1 

2fr 



de 



-Tr 



(20) 



where A(e) = i \G^ + '(e) — '(e)] is the total spectral function. It can be shown that the equilibrium torque 
always points out of the plane spanned by fi Ll and ri Ljv = , the magnetization directions of the two leads. This is 
expected, since the equilibrium contribution to the torques between the leads is a simply due to exchange coupling 
mediated by the paramagnetic system. The dynamics induced by this coupling conserves total energy and hence 
leaves il Ll • f2 LjVa invariant. Note that spin currents are even under time reserval, unlike charge currents, and can be 
present in equilibrium. 

The nonequilibrium (transport) contribution to the torques between the two leads, is given by 



ds 1 
~dt 



i r*+\e\v de 



ncq 



2h 



(27T) 



TrlhTl^i TG{ + l(e)-G[;)(e)T -i T KZ\™ '(e) - h^\e)r A\ tl {e) \ , (21 



:(-), 



and has both an out-of-plane and in-plane contribution. 
The former corresponds to a transport-modified electron- 
mediated exchange torque between the leads. The latter 
corresponds to the spin transfer torque between the fer- 
romagnetic leads. (This contribution was also calculated 
within the context of the Green's function formalism in 
Ref. [2l| for the case of two ferromagnetic leads.) Note 
that the in-plane spin transfer torque is present only 
at finite bias \e\V, which reflects the fact that energy- 
conservation violating torques can only be present in a 
nonequilibrium situation. Finally, note that this spin 
transfer torque acts on the leads and should not be con- 
fused with the torques that act on the local magnetization 
in the system that we calculate in the next section. 



III. APPLICATIONS 

The ferromagnetic voltage probe model discussed in 
the previous section is a flexible tool which can be used 
to model a wide variety of potentially interesting spin- 
tronic device geometries. We use it in this section to 
model both spin-independent and spin-dependent inelas- 
tic scattering. 

In the presence of an exchange field which defines 
a preferred direction in spin space, we can recognize 
(at least) three different length scales that are relevant 
to spintronic device functionality. The phase-coherence 
length is the length over which electronic quasiparti- 
cles maintain phase coherence. The spin-flip scattering 
length or spin diffusion length L a { is the average distance 
travelled along the channel between spin-flip scattering 
events. This quantity usually controls the length scale 



on which the magnititude of the magnetization along 
the exchange field direction recovers local equilibrium, 
and therefore also the paramagnetic spacer layer thick- 
ness scale at which giant magnetoresistance effects are 
strongly attenuated. The spin-coherence length L sc is the 
length scale over which components of the spin-density 
transverse to the exchange field can be maintained. The 
phase-coherence length is limited by spin-dependent in- 
elastic scattering processes. In addition to this inelastic 
contribution, the spin-coherence length in two and three 
dimensional ferromagnetic metals tends to be short be- 
cause of destructive interference due to the phase differ- 
ence of the spinor components in transverse conduction 
channels^ with different Fermi velocity components in 
the current flow direction. As we discuss below, phase 
coherence can have a strong influence on the properties 
of a spintronic circuit, especially in circuits containing 
antifcrromagnetic elements. 

The solution of transport Green's function models 
with ferromagnetic leads and non-collinear magnetiza- 
tion presents a number of numerical challenges, even 
for the one-dimensional case we consider here, because 
the spin and charge currents in the leads depend nonlin- 
early on both the magnetization direction and the mag- 
nitude of the exchange spin-splitting in each lead. When 
we wish to model a spin-independent inelastic scatterer, 
both the magnitude and the direction in each lead have 
to be carefully adjusted to achieve the zero-spin-current 
condition. Indeed the case of spin-independent inelas- 
tic scatterers can be very relevant to experiment since 
phonon-scattering can be dominant in experimental sys- 
tems. To model circumstances in which phonon scatter- 
ing is dominant (L s f 3> L^), we must deal with the com- 
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FIG. 3: Ferromagnetic GMR ratio as a function of the spacer 
thickness for various inelastic scattering lengths. 
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FIG. 4: Spin transfer torque for the case L B { — normalized 
to J/|e|, where / is the charge current, as a function of 0. The 
inset shows the conductance in units of \e\/(2nh). For this 
calculation the spacer thickness is taken equal to 5 sites. 



plications associated with preventing spin-relaxation in 
the voltage probes, a requirement that is especially trou- 
blesome for non-collinear magnetization configurations. 
Inelastic scattering off magnons, which does flip spins, 
can also be important however. The simplest case to 
examine numerically is one with paramagnetic, i.e. spin- 
isotropic, voltage probes. Calculations with this voltage 
probe model are relatively simple even for non-collinear 
magnetizations, since in this case only the charge current 
into the probes is set to zero. This limit corresponds to 
L s f = and can be interpreted as representing the case 
in which quasiparticle scattering from magnons is domi- 
nant. 

We start by considering the giant magnetoresistance 
(GMR) ratio which is easier to evaluate because we need 
to consider only parallel and antiparallel configurations 
which are both collinear. We define the GMR ratio as 



A. Ferromagnetic Metal Nanostructures 

We report results for a spin- valve ferromagnetic nanos- 
tructure, which consists of two ferromagnetic elements 
separated by a paramagnetic spacer (See Fig. This 
structure is known as a spin- valve because the current can 
be sharply reduced under favorable circumstances when 
the two ferromagnets have opposite orientations. In our 
toy model we choose A|/Js = 0.5 and independent of 
position, i.e. the site index j, in the ferromagnetic parts 
of the system. In the paramagnetic parts the exchange 
potential is set to zero. Both ferromagnets are 4 sites long 
in our tight-binding model. The other parameters used 
for the calculations reported on here were ef/Js = 1-8, 
Js = Jhj and V^f. 



Gp — Gap 



(22) 



1. GMR ratio for L s{ > L 4 



where Gp is the conductance for the parallel configura- 
tion and Gap for the antiparallel one. For antiferromag- 
nets, the designations parallel and antiparallel refer to 
the moment directions on the two sites adjacent to the 
paramagnetic spacer. 

In our model we define the phase coherence length as 
the product of the Fermi velocity and the inelastic scat- 
tering rate: 



4Jsa sin (hp a) 



E.-T^ep)/^ 



(23) 



where fcp is the Fermi wave length of the system. This 
definition is motivated the fact that the inelastic scatter- 
ing time is given roughly by 2/Tr J2j ^H £ f)/^b 



To model a finite phase coherence length due to spin 
independent inelastic scattering, the phonon scattering 
case, we choose L s f = oo. We then have to require 
that both spin and charge currents into the voltage 
probes vanish. As explained before, considering the 
non-collinear case for this situation is very difficult nu- 
merically. However, calculation of the GMR ratio re- 
quires comparing parallel and antiparallel configuration 
in which the spin currents with polarization out of the 
plane are zero. This observation makes calculating the 
GMR ratio numerically tractable. 

In Fig. 13 we present results for the GMR ratio as a 
function of spacer thickness for various inelastic scatter- 
ing lengths. The GMR ratio exhibits oscillations as a 
function of spacer thickness which have been observed 
experimentally- 6 . For the phase coherent case these os- 
cillations persist up to arbitrary large spacer thickness in 
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our one-dimensional model. For finite the oscillations 
are damped and the GMR ratio saturates to a nonzero 
value for a spacer of thickness much larger than L^. This 
is expected because for a ferromagnet the transmission 
and reflection probabilities are spin dependent, and not 
only the amplitudes. 

2. GMR ratio and spin transfer torques for = L s i 

For the magnon scattering case, = L s i, imple- 
mented by making the floating voltage probes param- 
agnetic and requiring the charge current to be zero, we 
are also able to consider the non-collinear case. The spin 
transfer torques are due to the net misalignment of elec- 
tron spins with the local magnetization in the transport 
steady stated, and calculated from 

r 8t = £^fif x( S ) x . (24) 

3 

Here, (s)^ denotes the electron spin component that 
points out of the plane spanned by the two magneti- 
zations. The sum over sites in the above expression is 
restricted to the ferromagnet for which we want to cal- 
culate the spin transfer torques, which in our case is by 
definition the right "downstream" ferromagnet. In Fig. 0] 
the spin transfer torque per current is shown for various 
as a function of the angle 9 between the two ferro- 
magnets. The inset of this figure shows the conductance 
as a function of this angle. Clearly, both GMR and spin 
transfer torques are suppressed for decreasing L$, as ex- 
pected in the L s f = case. 

B. Antiferromagnetic Metal Nanostructures 

In this subsection, we consider a system which con- 
sists of two antiferromagnets separated by a paramag- 
netic spacer. Hence, we take A| = A(-1) J 'in the mag- 
netic parts of the system, and zero in the paramagnetic 
parts. We take the parameters A/Js = 1. The other 
parameters used to obtain the results reported on here 
were e-p/Js — 1-8, Js — Jhj, and V^. The calculations 
presented in this subsection were performed on nanos- 
tructures in which each antiferromagnets has 30 sites. 

1. GMR ratio for L s{ > 

In Fig. [5] the antiferromagnetic GMR ratio is shown 
as a function of the thickness of the spacer, for various 
inelastic scattering length L$. Note that as the phase 
coherence length decreases, the oscillations in the GMR 
ratio with increasing thickness are suppressed. Unlike the 
case of a ferromagnetic spin valve, the GMR ratio goes to 
zero when the thickness of the spacer is much larger than 
the phase coherence length. This property demonstrates 
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FIG. 5: Antiferromagnetic GMR ratio as a function of the 
spacer thickness for various inelastic scattering lengths. 

that phase coherence is essential for the antiferromag- 
netic GMR effect. This observation is in agreement with 
qualitative arguments presented in previous workS. Note 
however, that the GMR decreases rather gradually with 
decreasing phase coherence length. 

2. GMR and current-induced torques for L$ — L s t 

One of the most important results in our previous work 
on antiferromagnets^ is that in the ballistic limit the out 
of plane spin density which is responsible for the current- 
induced torque is periodic with the period of the antifer- 
romagnet. For the tight-binding model used in this paper 
this property implies that this spin density component, 
present only in the transport steady state, is constant 
throughout the antiferromagnet. It is this property which 
makes current-induced torques very efficient in driving 
collective order parameter dynamics in uncompensated 
antiferromagnets in the absence of inelastic scattering. 

In Fig. we show the out of plane spin density for the 
angle 9 = it /2 between two moments in the antiferromag- 
nets opposite the spacer. This calculation was performed 
for a spacer thickness of 10 sites. In agreement with our 
previous results we find that for = 00 the out of plane 
spin density is constant. For nonzero the out of plane 
spin density decays away from the spacer into the anti- 
ferromagnets. 

In Fig. [7| we show the current-induced torque per cur- 
rent acting on the right antiferromagnet as a function of 
the angle 9 between the two moments facing the para- 
magnetic spacer. It is important to realize that a single 
domain antiferromagnet is driven by a staggered torque. 
Hence, in this case we have to calculate the staggered 
torque given by 

r st ^X>l) J 'Afifx( S > x . (25) 

3 

As expected, both GMR and current-induced torques are 
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FIG. 6: Spin density component out of the plane spanned by 
the moments of the antiferromagnets opposite the paramag- 
netic spacer layer for various inelastic scattering lengths. In 
this calculation = L E f. 
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reduced by inelastic scattering. It is somewhat surprising 
that the current-induced torques at = 20 are larger 
than at = oo. However, this is most likely because 
the current-induced torque is normalized to the current 
and because of the one-dimensional model under con- 
sideration. We checked that in the limit L$ — + the 
current-induced torques always vanish. 



IV. DISCUSSION AND CONCLUSIONS 

The most important conclusion of this paper is that the 
introduction of inelastic scattering does not immediately 
destroy magnetoresistive and current-induced torques in 
nanostructures with two antiferromagnetic elements sep- 
arated by a paramagnetic spacer. On the contrary, we 
find that the GMR ratio goes to zero smoothly for spacer 



thicknesses larger than the phase coherence or inelastic 
scattering length. (Ab-initio calculations for Cr/Au/Cr 
antiferromagnetic spin- valves^ demonstrate that there is 
an additional interface related contribution to antiferro- 
magnetic GMR which is not captured by our toy model 
and is not limited by inelastic scattering in the spacer 
layer.) Since typical inelastic scattering lengths at room 
temperature (~ lOnm) are much larger than the mini- 
mum paramagnetic spacer layer thicknesses required to 
make coupling between magnetic layers insignificant, in- 
clastic scattering does not impose any practical limitation 
on GMR and current-induced torques in thin film based 
antiferromagnetic metal nanostructures. 

Unlike the ferromagnetic case, current-induced torques 
in antiferromagnets do not follow simply from the ap- 
proximate conservation of total spin angular momen- 
tum. Toy model calculations, as well as ab initio calcula- 
tions for realistic systems^,, suggest that both GMR and 
current-induced torques tend to be somewhat weaker for 
very thin films in the antiferromagnetic case compared 
to GMR and spin transfer torques in ferromagnets. It 
seems natural to associate this property with the absence 
of grounding in a simple robust conservation law. A sur- 
prising and interesting property of ST physics in antifer- 
romagnets is the property that torques act throughout 
the entire volume of the antiferromagnet in the phase 
coherent case. This property helps to compensate for 
the tendency toward somewhat weaker effects in very 
thin films. When inelastic scattering is introduced, spin- 
torques again act only over a finite thickness ~ within 
the antiferromagnets. Still the length scale over which 
the torque acts should be ~ lOnm, much longer than 
the atomic length scale over which spin torques act in 
ferromagnetic spin-torque structures. 

As explained in earlier work^, ST effects in antifer- 
romagnets are expected to be strongest when moment 
directions alternate in the direction of current flow. In 
the one-dimensional single channel model studied here, 
this is the only possibility. For real three-dimensional 
antiferromagnets this property requires particular orien- 
tations of current with respect to the crystal axes, or in 
thin film structures particular growth directions for the 
antiferromagnetic films which yield alternating ferromag- 
netic layers. It seems likely that the greatest obstacle to 
achieving interesting GMR and current-induced torques 
in purely antiferromagnetic circuits may be the difficulty 
in growing antiferromagnetic thin films with individual 
ferromagnetic layers that are not strongly compensated. 
Studies of the influence of magnetic disorder within fer- 
romagnetic layers on GMR and ST physics require multi- 
channel Green's function calculations and will be the sub- 
ject of future work. 
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